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A description of mesoscopic fluctuations of the pairing gap in finite-sized quantum systems based 
on periodic orbit theory is presented. The size of the fluctuations are found to depend on quite 
general properties. We distinguish between systems where corresponding classical motion is regular 
or chaotic, and describe in detail fluctuations of the BCS gap as a function of the size of the system. 
The theory is applied to different mesoscopic systems: atomic nuclei, metallic grains, and ultracold 
fermionic gases. We also present a detailed description of pairing gap variation with particle number 
. . . for nuclei based on a deformed cavity potential. 
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C3 ' INTRODUCTION 

hj : 

Bohr, Mottelson and Pines were first to apply the Bardeen, Cooper, Schrieffer (BCS) theory of superconductivity 
to finite size systems, namely to describe pairing in atomic nuclei [1[. A consequence of the finite size of the system 
i— i, is the appearance of shell structure. This implies fluctuations of the pairing gap as a parameter, like the particle 
i^h ■ number, is varied. In this contribution we shall discuss how these fluctuations can be described in a semiclassical 
I , theory. 

We first discuss pairing in nuclei obtained from the odd-even mass difference. Pairing gaps calculated from different 
mass models are compared, both with respect to average gaps and to fluctuations. In the next section a semiclassical 
theory of fluctuations of the BCS pairing gap Q is presented. This provides analytic expressions of fluctuations of 
the pairing gap, where the dynamics of the underlying classical system (chaos/order) is an important parameter. The 

t-H ■ theory is first applied to pairing fluctuations in nuclei, where a detailed comparison with data is performed. We also 
utilize explicit periodic orbits, taken from a deformed cavity model, to describe the detailed variation of the pairing 
gap with particle number. Nano-sizcd metallic grains are also studied, where, due to a lack of experimental data, we 
compare our theoretical results to other existing numerical results. Interaction strength, external potential, as well as 

(T) ■ the number of particles can be experimentally tuned in ultracold fermionic gases, and we discuss the size of pairing 
* , fluctuations for such systems, as obtained from our theory. 

o : 

oo . 

O ODD-EVEN MASS DIFFERENCE IN NUCLEI 

>' 

The systematic difference between the ground-state mass of odd and even nuclei constitutes an important indicator 
r> ' of pairing in nuclei. The pairing gap can be calculated from binding energies, B, utilizing the three-point measure 

A 3 (M) = B(M) - ~[B{M + 1) + B(M ~ 1)], (1) 

where M is the neutron N or proton Z number. In the presence of other possible interactions, this quantity has been 
shown to be a good measure of pairing correlations Q provided M is taken as an odd number. In that case, it is 
easy to see that there is no contribution from the mean field in A3, while if M is taken as an even number, an extra 
contribution to A3 from the single-particle levels appears in the extreme single-particle model as ^(e^+i — e^), where 
ei is the last occupied single-particle level. In Fig. Q] we see the systematic difference in the average of A 3 when M is 
an odd and even number. Restricting to A 3 -values obtained from odd numbers of M , a fit of the pairing gap gives 

A = MeV - (2) 

where A = N + Z is the total number of nucleons. If, however, also cases with M is an even number are included, 
the usually employed pairing gap value, 

A = ^2 MeV > (3) 
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FIG. 1: Average pairing gap, A, versus particle number A, as obtained from Eq. ((TJ using measured masses [j. The averaging is 
performed in a region of fixed mass number, A, over all available isotopes. The particle number M (neutron or proton number) 
is either odd (lower black dots), even (upper red dots), or both even and odd (middle blue dots). The curves, 2.7 /A 1/4 MeV 
(lower line) and 12/ A 1 ^ 2 MeV (upper line), are obtained from fits to cases with odd M, and both odd and even M- values, 
respectively. 
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FIG. 2: Average pairing gap versus mass number A using the three-point measure, Eq. Q, with odd M. Binding energies arc 
obtained from different theoretical mass models, Moller et al [f| (black dots), Samyn et al (red triangles pointing up) [(J and 
Duflo et al (magenta triangles pointing down) [?J . The theoretical curves are compared to experimental results (black dots) , cf 

Fig. m 



gives a better fit. The fitted difference A3(even M)-A3(odd M) = ^(e^+i — e*) is found to be inversely proportional 
to the mass number, A, and vary as 50/A MeV. The pairing gap as defined by Eq. (JTJ) can also be extracted 
from theoretical mass models, and in Fig. [2] we compare the average value of A3 versus particle number for three 
different mass models. Two of them are based on mean field theory. The first one is a non-self-consistent macroscopic- 
microscopic model 5|, the second one is a self-consistent calculation based on Skyrme-Hartree-Fock-Bogoliubov @ 
while the third one 7j is a shell-model based calculation with parameterized monopole and multipole terms. Also the 
experimental mean value of the pairing gap is shown in Fig. [2l and it is seen that all three mass models give similar 
results in good agreement with experimental numbers. If, however, not average values, but pairing gaps obtained from 
all nuclei are shown (Fig. [3]), it is clear that the average values give a poor description of the result. The fluctuations 
(RMS value) of the pairing gaps are indeed very large, and exhibit different variation with particle number for the 
different mass models, see Fig.fU All mass formula show the same tendency of decreasing fluctuations with increasing 
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FIG. 3: Pairing gaps extracted from Eq. ([TJ with odd M versus mass number A from experimental data, and three different 
theoretical mass models. Each black point corresponds to a nucleus. The solid red lines show the average pairing gap as shown 
in Fig. H 
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FIG. 4: Root mean square (RMS) values of pairing gaps A obtained from Eq. ([T]) using odd M-values for three theoretical 
mass models, and experimental data. See caption of Fig. [2] for notations. 



mass number, as is also seen from experimental masses. However, considerably larger pairing fluctuations are seen in 
the mass model by Samyn et al [6(, where, particularly for large mass numbers, almost three times larger fluctuations 
are obtained, as compared to pairing gap fluctuations obtained from measured masses. The mass model based on 
the shell-model gives systematically too small fluctuations, while the mass model by Moller et al [f| gives pairing 
fluctuations closest to experimental data. 

It is thus clear that the fluctuations of the pairing gap is a most important property, and in the next section we 
shall present a semiclassical theory for pairing fluctuations based on periodic orbits [2| . 
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PERIODIC ORBIT DESCRIPTION OF PAIRING FLUCTUATIONS 



The many-body Hamiltonian, 

H = Hi + Vp a i r = ^ eka[a k - G^^a\a\aia t} (4) 



k( 



incorporates a one-body part, Hi (typically obtained from a deformed mean-field), and a two-body pairing interaction 
between time-reversed states, V^air- All two-body matrix elements are assumed to take the same value, G (seniority 
interaction). In the mean-held approximation in pairing space a pairing gap, or pairing "deformation", 



is determined by the BCS gap equation 



that can be written as 



2 



A=(g5>MV (5) 



(0) 



G = \ ^(e k -\f-^ 



2 f L p{e)de 



G J_ L Ve 2 + A 2 ' 



(7) 



where p(e) is the single-particle level density, and we have put the Fermi energy, A, to zero. The energy cut off is at 
±L. 

Following semiclassical approaches, the pairing gap as well as the single-particle density of states are divided in a 
smooth part and a fluctuating part, A = A + A and p = p + p, respectively. In the weak coupling limit A«L, the 
smooth part of the gap is given by the well known solution A = 2Lexp(— 1/pG) (see e.g. Ref. [8J). In semiclassical 
theory the fluctuating part of the density p can be calculated from purely classical properties [j| , 

oo 

p(e) = 2 A p>r cos(rS p /h + v Pir ), (8) 

p r—l 

where the sum is over all primitive periodic orbits p (and their repetitions r) of the classical underlying effective single- 
particle Hamiltonian, H\. Each orbit is characterized by its action S p , stability amplitude A p>r , period r p = dS p /de 
and Maslov index v p ^ r (all evaluated at energy e). Assuming Ac A and A<L gives after some algebra Q 

~ A °° ( rS \ 

A = 2-Y J J2 A P> rK °( rr p/ T ^ C0S [ ^h +u P> r ) ' ( 9 ) 

P p r—l ^ ' 

where all classical quantities involved are evaluated at the Fermi energy. Kq(x) is the modified Bessel function of 
second kind, and 

(10) 



2ttA 



is a time corresponding to the pairing gap, that we may call the pairing time. Since Kq(x) oc exp(— x)/^/x for 
x ^> 1, the Bessel function exponentially suppresses all contributions for times r ^> r A (making the sum in Eq. ([9|) 
convergent) . 

Since the value of the actions depend on the shape of the mean-field potential, Eq. ^ predicts generically fluctu- 
ations of the pairing gap as one varies, for instance, the particle number, or the shape of the system at fixed particle 
number. The fluctuations result from the interference between the different oscillatory terms that contribute to A. 
When the motion is regular (integrable) , continuous families of periodic orbits having the same action, amplitude, 
etc, exist. The coherent contribution to the sum ([9]) of these families of periodic orbits produces large fluctuations. 
In contrast, in the absence of regularity or symmetries, incoherent contributions of smaller amplitude coming from 
isolated unstable orbits are expected for chaotic dynamics. 
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FIG. 5: The two functions Fa and Fi, see Eq. (|14p. versus the argument D. 



The second moment of the fluctuations may be obtained from Eq. <j9j) as 



A 



r 2 , 

'e jo 



A > =V / dr^(r/r A )lf(r), (11) 



where t h = h/ 5 is the Heisenberg time ((5 = is the single-particle mean level spacing at Fermi energy), and K(t) 
is the spectral form factor, i.e. the Fourier transform of the two-point density-density correlation function [10| . 

The structure of the form factor K(t) is characterized by two different time scales. The first one, the smallest 
of the system, is the period T m j n of the shortest periodic orbit. The form factor is zero for r < r m j n , and displays 
non-universal (system dependent) features at times r m ; n ~ r <C t h . As r further increases, the function becomes 
universal, depending only on the regular or chaotic nature of the dynamics, and finally tends to t h when r ^> t h . The 
result of the integral (jlip thus depends on the nature of the dynamics, and on the relative value of r A with respect 
to T m i„ and t h . In the simplest approximation, all the non-universal system-specific features are taken into account 
only through T m - ln [HI], and one can write K(t) — for r < r m - m and, for r > r mul , K(t) = t h for integrable systems 
and K(t) — 2r for chaotic systems with time reversal symmetry. We assume a generic regular system; the analysis 
does not apply to the harmonic oscillator, whose form factor is pathological. 

This finally gives the expressions for fluctuations of the pairing gap (normalized to the single-particle mean level 



spacing), a = ^/(A 2 )/^, assuming regular dynamics 

7T A 



and assuming chaotic dynamics, 



where we have introduced the function 0] 



The argument D is denned as 



ai s = ~jF (D) , (12) 



^ 2 h = (D) , (13) 



F n (D) = l- J ° ° . (14) 
J Q x n K^{x)ax 



r A g 6 



where the parameter g is the ratio between the Heisenberg time and the time of the shortest periodic orbit, 

th 



(16) 
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FIG. 6: Fluctuations of the pairing gap as a function of the mean value for mesoscopic systems (log-log scale; all quantities 
normalized with the single-particle mean spacing). Regular and chaotic dynamics are shown by blue and red lines, respectively. 
The dashed curves correspond to different values of the dimensionless conductance, g, and the limiting case of g — > oo is shown 
by solid lines. The results are valid in the pairing regime A/5 > 1. Applications to Nuclei and Nano-grains are marked out. 
Ultracold atomic gases may be controlled to appear in major parts of the figure. From Ref. Q. 



This parameter is often called "dimensionless conductance". It expresses the energy range, L max (= g) (in dimen- 
sionless units; energies are divided by the mean energy spacing, $), in the single-particle spectrum where spectrum 
fluctuations show universal properties [lfl ]. For a system where the corresponding classical dynamics is regular/chaotic 
(with time reversal symmetry), the statistical properties of the one-body energies, of Eq. (H|) are described by Pois- 
son/GOE (Gaussian Orthogonal Ensemble) statistics, respectively. The limit g — ► oo (i.e. D—0), corresponds to a 
universal situation when the full spectrum corresponds to pure GOE (if chaotic) or Poisson (if regular) statistics. 

The dimensionless parameter D can also be expressed as the system size, 2R, divided by the coherence length of 
the Cooper pair, £ = hv F /(2A), where v F is Fermi velocity, 

D = 2R/&. (17) 

The Cooper pairs can thus be considered as restricted by the system size if D < 1. 

In Fig. [5]we show the two functions Fq and F\ of Eq. lfT4|) versus the argument D. For small values of D, Fo ss 1 
and F\ w 1, and a universal behavior appears, i.e. the fluctuations do not depend on the system properties, only the 
dynamics, implying Poisson and GOE statistics for the full spectrum in the regular and chaotic cases, respectively. 
In the other limit when D becomes large, Fq and F\ go to zero and all pairing fluctuations disappear. This is, for 
example, the situation in bulk systems with a large number of particles. 

In Fig. [6] pairing gap fluctuations are shown versus the pairing gap for different values of the dimensionless conduc- 
tance, g, for regular as well as for chaotic dynamics. The plot covers a large range of parameter values and is shown 
in log- log scale. Pairing gaps of the order of the mean level spacing or smaller, A < 5, (ultrasmall regime), are not 
treated by the present theory and corresponding region is shaded in the figure. In this region the Anderson condition 
[r3 | implies no BCS pairing. As mentioned above, non-universal behavior (deviating from GOE or Poisson) appears 
when D > 0, corresponding to finite values of the dimensionless conductance, g. 
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FIG. 7: Fluctuations of pairing gap versus average pairing gap for different values of D. Regular dynamics is assumed. D = 0.3 
approximately corresponds to nuclei. Fluctuations and pairing gaps are expressed in units of the mean level spacing, 8. 



FLUCTUATIONS OF PAIRING GAP IN NUCLEI 



In a previous section we studied pairing gap fluctuations in nuclei as obtained from the odd-even mass difference, 
where binding energies were taken from different mass models as well as from data, see Fig. |4] We now would like to 
use the semiclassical theory developed in previous section to compare pairing gap fluctuations, and compare to these 
results. Notice that the theory does not contain any parameters. Once the system is defined, and the variable D 
has been determined, the fluctuations only depend on the corresponding classical dynamics, being regular or chaotic. 
In nuclei ground states the dominating dynamics is expected to be regular, although elements from chaos may be 
present, see Ref . [l3l| . 

To evaluate the root mean square value (RMS) of the pairing fluctuations from the theoretical expressions, Eqs. (fT2"|) 
and (|13[) . we only have to determine the parameter D. The size of the nucleus, 2R = 2x fm, and size of the 

pairing correlation length, £ = hv F /(2A) = 11.3A- 1 / 4 fm, give (using Eqs. © and (JUJ)) D = 2i?/£o = S-22A 1 / 12 = 
0.27 — 0.33 for mass numbers in the interval, ^4=25-250. The dependence of the pairing fluctuations, assuming regular 
dynamics, with the (average) pairing gap is shown in Fig. [7] for the three cases, D—0, 0.3 and 1. If the pairing 
correlation length is smaller than the system size {D > 1) pairing fluctuations are quite small. Largest fluctuations 
appear in the universal limit when D=0. The small values of D for nuclei (D sa 0.3) implies that the Cooper pairs are 
non-localized. The pairing gap fluctuations are thus substantial, and are about half the value at the universal limit. 

We may compare the fluctuations to the experimental pairing gap fluctuations by inserting the above values for D, 
and the average values of A from Eq. ([2]), in Eqs. (| L 2[) and (fT3)) . assuming regular or chaotic dynamics. The resulting 
curves are compared to the experimental one in Fig. [5] Note the good agreement between the theoretical pairing 
gap fluctuations assuming regular dynamics, and the experimental curve, both in the overall amplitude and in the 
A-dependence. In Ref. [l3j it was discussed the possibility that the dynamics of the nuclear ground state is mixed 
regular and chaotic. Making this assumption in the calculation of fluctuations of the pairing gap results in a curve 
that is very close to the purely regular curve in Fig. [8] That is, fluctuations of the pairing gap cannot distinguish if 
there is a chaotic component in the nuclear ground state. 



SHELL STRUCTURE IN PAIRING GAP FROM PERIODIC ORBIT THEORY 



Shell effects in nuclei are also seen in the pairing gap. One may go beyond a statistical description, and use Eq. ([9]) 
to obtain a detailed description of the variation with neutron or proton numbers. For that purpose, we assume for 
the nuclear mean field a simple hard-wall cavity potential. The shape of the cavity at a given number of nucleons is 



fixed by minimization of the energy against quadrupole, octupole and hexadecapole deformations To simplify, 



we take N = Z. The periodic orbits of the spherical cavity (a few short orbits are shown in Fig. [9]) are used in Eq. 
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FIG. 8: RMS values of pairing gap fluctuations in nuclei versus mass number A, obtained from measured masses (dashed line 
connected by filled circles). Semiclassical calculations of the fluctuations assuming chaotic (Eq. ((13}) and regular (Eq. (112[) ) 
dynamics are shown by red dashed and blue solid lines, respectively. From Ref. @] 
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FIG. 9: A sample of classical periodic orbits of the cavity potential. The orbits shown correspond to the index (v,w) — (2,l), 
(3,1), (4,1), (5,2), (7,2) and (7,3), counted from upper left to lower right figure. 

with modulations factors that take into account deformations and inelastic scattering [l5| . This gives 

A = — — A vw M vw (x) K x (£ vw ) Kg \ J,- W ^ ) sin(k F R£ vw +v vw ^-), (18) 

where M vw (x) is a modulation factor for perturbative deformations, n x (( V u) a modulation factor for inelastic scat- 
tering, and x stands for the three considered deformation degrees of freedom, quadrupole (£2), octupole (e 3 ) and 
hexadecapole (£4). The summation is carried out over the two indices (v,w) including the periodic orbits shown in 
Fig.[5J We set the average of A to zero, as was done with the experimental data. In Fig. [TO] we compare the theoretical 
result A(iV) to the experimental value averaged over the different isotopes at a given N. The agreement is excellent; 
the theory describes all the main features observed in the experimental curve. 



PAIRING FLUCTUATIONS IN NANO-SIZED METALLIC GRAINS 

Experiments in the 90's have explored the superconducting properties of nanometer scale aluminum grains [lij ]. 
Irregular shape of the grains implies chaotic dynamics, and energy levels are found to follow GOE (there are no further 



symmetries than time- reversal), see Ref. 17] . The existence of a superconducting gap was demonstrated in the regime 
A > S, whereas no gap was observed when A < 8. The transition occurs around N ~ 5000, where N is the number of 
conduction electrons in the grain. The N dependence of the average gap A is poorly understood. We will adopt for 



grains the thin-film value A w 0.38 x 10" 3 eV jl7|. The mean level spacing is 6 = (2E F )/(3N) w 2.1 /AT eV, whereas 
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FIG. 10: Nuclear pairing gaps for neutrons. Experimental isotope sequences are connected by thin solid lines. The blue dashed 
line shows average pairing gaps. Calculations from the cavity model are shown by the solid red line. Average behavior (Eq. ([2J) 
is shown by solid black line. Data from Ref. Q. 



the dimensionless conductance, g p=> 2.6-/V 2 / 3 . Eq. (fT5|) gives D fa 4.4 x lO^N 1 / 3 , which ranges from 0.05 to 0.02 
when N varies between 10 3 and 10 5 . This means that the variance will be close to its "universal" value (GOE limit) 
obtained by setting Fi(D) = 1 in Eq. fl3|). namely [3 
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1/27T 2 



(19) 



The typical range of variation of pairing gap fluctuations for nano-grains is marked out by a grey area in Fig. [6l 
There are no explicitly measured pairing gap fluctuations for nanosized metallic grains to compare the present 

theory. We may, however, compare to other independent calculations, namely the condensation energy calculated of 

chaotic grains [20j , and variation of the pairing gap with particle number in a regular cubic shaped grain 21 1 . 

The BCS condensation energy is defined as the total energy difference between the paired and the unpaired system. 

With our choice of gap equation |J7J| it is written as 



E c =E tot (A)-E tot (A = 0) = J p{e)ev 2 {e)de 



A 2 
4G 



p{e)edx 



(20) 



where v 2 (e) = ^ ^1 — ^ J ■ Inserting the semiclassical approximations of p = p+p and A = A + A and expanding 
Eq = Eq + Ec to lowest order in fluctuating properties, assuming Eq <§C Eq , gives 



E c = -2 ^2 ^2 A p,r [Qp, r + qK (rT p /T & )} COS 

p r—1 



v. 



V ■ r 



(21) 



where Kq is the modified Bcssel function of second kind, q = A 2 ( 1 — jjio) a '"' 



Q P ,r = J cos (^) e ^1 



Ve 2 + A 2 



de 



(22) 
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FIG. 11: Fluctuations of the condensation energy as a function of the mean level spacing of the system. The blue solid and 
red dashed lines show the numerical calculations of regular and chaotic fluctuations according to Eq. ((23}. The squares show 
the calculations of ref. 12011. 



The second moment, er 2 = (E c 



is thus obtained as 



2 

k 1 



[Q(r)+qK (r)} 2 K(r)dT, 



(23) 



where K(t) is the spectral form factor. Using in addition the estimates G — 0,224 and L = A e V - 224 [2(| we show 
in Fig. [TT] the fluctuations of the condensation energy, ac, versus the level distance in the system, assuming regular 
and chaotic motion. As mentioned, there is no experimental data to compare to, but we may compare to a numerical 
calculation of the ac by Sierra et. al. [20J where they use the Richardson's solution of the pairing problem and random 
matrix theory (GOE) for generating the spectrum. We see that our semiclassical theory, assuming chaotic dynamics, 
agrees well with the random matrix calculation of Ref. [20| ■ The random matrix model applied in [201 ] happens to be 
a reasonable approximation for the considered nano-grains, since D ss 0, i.e. the nano-grains have properties which 
are close to the universal limit. 

Garcia-Garcia et al [2l| studied shell effects on the pairing gap in a cubic geometry as a function of the mean 
pairing gap, A/ 5. The cubic geometry implies classically regular motion, and we may thus use Eq. (|12p above to 
calculate the fluctuation of the paring gap, normalized to the mean pairing gap A instead of the mean level spacing 

<5, a rcg S/ A = \/f 4) where we have set F = 1 since for small grains D w 0. This gives a very good agreement to the 



pairing fluctuations calculated in Ref. [2l|. It is interesting to note that by making the system chaotic, the pairing 
fluctuations decreases substantially and become (see Eq. (TP3"]Q , a^S/A 
A/6 = 20. 
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PAIRING FLUCTUATIONS IN ULTRACOLD FERMIONIC GASES 



Recently, a large interest has emerged in studying trapped atomic gases of bosons and fermions. The gases are 
ultracold and dilute, and provide the possibility to study new quantum phenomena in the physics of finite many- 
body systems. The number of neutrons of the confined atoms determines the quantum statistics; odd number implies 
Fermions, and even number implies bosons. Studies of Bose condensates can be done for the bosonic gases, and studies 
of quantum phenomena including superfluidity and the transition to a Bose-Einstein condensate can be conducted 
for the Fermi gases. 

Since the atom-atom interaction is short ranged and much smaller than the interparticle distance, the atom-atom 
interaction can be approximated by the (5-interaction, 



V(n - ra) 



4^5(3) ( ri 



ra), 



(24) 
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where a is the s-wave scattering length, that can be externally controlled in size and even in sign through the Feshbach 
resonance. Also the confinement potential can be externally controlled to create regular as well as chaotic dynamics. 
Ultracold fermionic gases thus provide excellent conditions for theoretical as well as experimental studies of pairing 
properties, including pairing fluctuations. Since both particle number and interaction strength are experimentally 
controlled parameters, the fluctuations may appear in major parts of Fig. [SJ 

We estimate 5 = (2E F )/(3N) and g = ±(3A0 2/3 ; in the dilute BCS region A = (2/e) 7 / 3 E F exp (-iv/2k F \a\) [13], 
with k F the Fermi wavevector, giving 

D = ^ =2^(2/e) 7/3 (3A0 1/3 exp(-7r/2/Ma|). (25) 

Recent experiments using Li 6 reach k F \a\ — 0.8 [23T |. implying negligible fluctuations for typical values of N ~ 10 6 . 
Reducing to k F \a\ — 0.2 and N — 10 4 yields for generic regular systems fluctuations that are on the same magnitude 
as the mean pairing gap, cr reg w 0.5A/5. 

SUMMARY 

In summary, we have presented a semiclassical theory that provides a generic description of fluctuations and shell 
structure of the BCS pairing gap in finite Fermi systems. These mesoscopic systems are generically dominated by 
system specific features not included in purely statistical models like GOE. Different possible regimes, as well as the 
influence of order/chaos dynamics, were investigated, in particular for the typical size of the fluctuations (Fig. 
The present theory provides analytic predictions, valid for a wide range of physical situations. It also compares very 
favorably with available experimental data. 
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